﻿ Filtre numerice transversale 1 Scopul lucrării Se studiază în domeniile timp şi frecvenţă, filtrele transversale în timp discret Prin particularizarea structurii acestor filtre se obţine mediatorul numeric alunecător Se studiază acest tip de sistem din perspectiva filtrelor adaptate 2 Filtre numerice transversale în timp discret Ecuaţia cu diferenţe finite care caracterizează filtrul numeric, care la excitaţia [ ]sn, ră ispunde cu semnalul []sn este: 0 NM [] []asnk bsnk−=− (1) 0ki k∑∑ 00kk== Condiţiile iniţiale se consideră nule Dacă singurul coeficient a nenul este a atunci filtrul implementat este un sistem nerecursiv Aceste k0 filtre au proprietatea ca răspunsul lor la impuls este de durată finită şi de aceea se mai numesc filtre cu ia care descrie un astfel de sistem este: răspuns finit la impuls (FIR) Pentru 1a= 0, ecuaţ M [] [ ]sn bsn k= (2) ∑ 0ki− 0k= Forma canonică I de implementare a sistemului descris de ecuaţia (2) este prezentată în figura 1 []sn i D D D b b 12 bb M0 + 0[]sn Fig 1 Forma canonică I de implementare pentru un sistem FIR Pe baza figurii se constată ca operaţiile de întârziere, cele de înmulţire cu constante şi cele de însumare se realizează pe direcţii transversale, de unde vine şi denumirea de filtre care se implementează în acest mod 1Relaţia (2) conduce la următoarea expresie a răspunsului la impuls al filtrului transversal: M − (3) [] [ ]hn b n kδ ∑ k= 0k= din care se poate deduce expresia răspunsului în frecvenţă: M Ω jk− ()HbeΩ= (4) k∑ 0k= 3 Mediatorul numeric de semnal Dacă coeficienţii b au valoarea 1M pentru 0, 1kM= şi în rest sunt nuli se obţine ecuaţia cu k− diferenţe finite: 1M− (5) [](1/)[]sn M sn k= ∑ 0i− 0k= Având în vedere că semnalul []sn se obţ 0ine din semnalul [ ]isn, mediind aritmetic ultimele M eşantioane ale acestuia, sistemul obţinut se numeşte mediator numeric Răspunsul său la impuls este: 1M− − (6) [] (1/ ) [ ]hn M n kδ ∑ 0= 0k= iar răspunsul său în frecvenţă este: 1M− (1)/2jk j M−Ω −Ω − = Ω Ω (7) ( ) (1 / ) (1 / ) sin( / 2) / sin( / 2)HMeMe MΩ= ∑ 0k= răspunsul la impuls şi modulul răspunsului în frecvenţă ale mediatorului numeric sunt Pentru 4M= prezentate în figura 2 h[n] 1/4 n -1 0 1 2 3 4 Fig 2 Răspunsul la impuls şi modulul răspunsului în frecvenţă ale mediatorului numeric, M=4 2Se constată că mediatorul este un filtru trece jos, care la anumite frecvenţe introduce atenuări infinite Poziţia pe axa frecvenţelor a acestora precum şi atenuarea minimă în banda de blocare pot fi fixate cu ajutorul constantei M Analizând ecuaţia (5), se constată că, sincron cu derularea temporală a eşantioanelor semnalului []sn, 0 are loc “alunecarea temporală” a mediei aritmetice a ultimelor M eş antioane ale semnalului [ ]sn De i aceea sistemul prezentat este numit mediator alunecător 4 Filtre numerice adaptive la semnale în timp discret Se consideră că trebuie transmisă secvenţa utilă []sn, de durata finită M, acoperită aditiv de zgomot i []nn În scopul reducerii zgomotului din semnalul: i [] [] []xnsnnn= (8) ii+ hn , acest semnal este prelucrat de către sistemul nerecursiv în timp discret cu răspunsul la impuls [ ] obţinându-se semnalul [ yn: ] (9) [] [] []yn s n n n= 00+ unde []sn reprezintă a semnalului [xn iar []nn componenta sa aleatoare ] 0 componenta utilă0 Problema este determinarea răspunsului la impuls [ hn al acelui sistem, care asigură la ieşirea sa un ] raport semnal pe zgomot maxim, la un moment de timp n Sistemul care satisface aceasta condiţ 0ie se numeşte filtru adaptat la semnalul [ ]sn Semnalele [ ]sn ş finită deci de putere iii 0[]sn sunt de durată nulă Raportul semnal pe zgomot la ieşire trebuie să depindă de timp, pentru a putea fi maxim la momentul n Din aceste motive formula de calcul al raportului semnal pe zgomot la ieş 0ire este: 2 []/ []SNR s n P n= 00 0n (10) unde cu []Pn s-a notat puterea zgomotului de la ieş 0nire Dar: 1M− [] []*[] [ ] [ ]s n s n hn s p hn p== − (11) ∑ 0ii⋅ 0p= Folosind inegalitatea Cauchy-Bunyakovsky-Schwarz, se poate scrie: 11MM−− ⎛ ⎞ 222⎛⎞ (12) [] [ ] [ ]sn s p hn p≤− ⎜ ⎟ ∑∑ 00i⎜⎟ 00pp== ⎝⎠⎝ ⎠ Egalitatea are loc în ultima relaţie dacă: 0, 1pM∀=− (13) [][]hn p as p− 0i= 3 Deci sistemul care răspunde la semnalul [ ]sn cu un semnal []sn, maxim la unde a este o constantă i0 momentul n, este acela al că spuns la impuls este: 0rui ră − (14) [] [ ]hn as n n=− 0i 00(1),nn M n=− dacă 1nM=− 0, atunci: − 0, 1nM= (15) [] [ 1 ]hn as M n=− i− Se consideră ca []nn este un semnal aleator de tip zgomot alb, ale că antioane nu sunt corelate cu 0rui eş , din toate sistemele cele ale semnalului [ ]sn ş ii sunt independente între ele Se poate demonstra că nerecursive cu răspunsul la impuls de durata M, acela care asigură la ieşirea sa un raport pe zgomot ia (14) 0SNR maxim este cel dat de relaţ De aceea, acest filtru adaptat la semnalul [ ]sn în ipoteza ca acesta este perturbat aditiv de zgomot alb i Dacă semnalul util de la intrarea filtrului adaptat are expresia: 1, 0, 1nM⎧=− ⎪ (16) []sn= i⎨ 0,rest⎪ ⎩ şi dacă constanta a este aleasa de valoare 1M, atunci expresia ră spunsului la impuls al filtrului adaptat la acest semnal este: 1/ , 0, 1Mn M⎧=− ⎪ []hn= (17) ⎨ 0,rest⎪ ⎩ Deci mediatorul alunecător este un filtru adaptat pentru semnalul definit de relaţia (16) Acest semnal este de durata finită şi amplitudine constantă Mediatorul numeric în timp discret poate fi folosit pentru îmbunătăţirea raportului semnal pe zgomot al semnalelor periodice în timp continuu 5 Desfăşurarea lucrării Obiectul acestei lucrări de laborator este un sistem cu schema bloc din figura 1, simulat in Matlab Semnalul util de la intrarea sistemului este o sinusoidă în timp discret, cu N eşantioane pe o perioadă (maxim 512) Numărul M, putere a lui 2, are valorile 2, 4, 8, 16, 32, 64 5 1 Folosind programul ”ftrans”, se vizualizează şi se reprezintă răspunsul la impuls al sistemului şi i [ pentru M de valori 2, 4, 8, 16, 32 şi 64 semnalele de intrare/ieşire[ ]] isn ş0sn 5 2 Se determină caracteristicile de frecvenţă ale sistemelor simulate la punctul 5 1 şi se reprezintă pe hârtie milimetrică Se compară cu rezultatul obţinut în Matlab 45 3 Se amestecă semnalul [sn cu un zgomot alb de putere cunoscută, obţinându-se semnalul [xn ] i] Raportul semnal pe zgomot la intrare SNRi se calculează împărţind pătratul valorii efective a semnalului [ ] isn la puterea zgomotului Se alege lungimea filtrului M=2, 4, 8, 16, 32, 64 Se calculează raportul semnal pe zgomot la ieşire i zgomot Se calculează îmbunătăţirea SNRo, împărţind pătratul valorii efective ale semnalelor [ ] 0sn ş raportului semnal pe zgomot Γ folosind raportul SNRo/SNRi Se fac operaţiile pentru celelalte 5 valori pe hârtie milimetrică dependenţa (MΓ ale lui M Se reprezintă ) 6 Exerciţii în MATLAB 6 1 Fenomenul Gibbs 6 1 1 Semnalul dreptunghiular poate fi reprezentat ca o serie de semnale sinusoidale sau cosinusoidale armonice Considerăm suma parţială a următoarelor semnale periodice: ∞ 411 41⎡⎤ , unde n este un număr impar ωω ω… =+++= () ∑ 000 0sin sin 3 sin 5 sinxtttt ntω ⎢⎥ ⎣⎦ = 135nnππ Exemplul următor de comenzi MATLAB foloseşte o buclă pentru a genera această sumă pentru orice număr impar: n=input('Introduceti un numar impar: '); omega 0=0:0 01:2*pi; x=0; for k=1:2:n; x=x+1/k*sin(k*omega 0); end x=4/pi*x; plot(omega 0,x), xlabel('omega 0') text(3 5, 7,['Suma a ', num2str((n+1)/2), ' unde sinusoidale']) Editaţi şi rulaţi secvenţa de mai sus pentru diferite numere impare 6 2 Transformata Fourier în timp discret 6 2 1 Să se calculeze transformata Fourier discretă a secvenţei: [8xn n nσσ=−−, pentru ][][ ] 1 020n≤≤ Să se reprezinte grafic partea reală, partea imaginară, modulul şi faza transformatelor Fourier discrete calculate n=0:20; x1=[ones(1,8), zeros(1,13)]; stem(x1) X=fft(x1,512); %vezi help fft Calculul transformatei Fourier discrete s-a efecuat într-un număr de 512 de puncte plot(abs(X)),grid; %vezi help abs Această reprezentare corespunde intervalului de frecvenţă [0,2π) Prin utilizarea comenzii fftshift se reprezintă în intervalul [-π, π) În acest scop se generează un vector cu pas linear care conţine în 5intervalul respectiv un număr de elemente egal cu lungimea transformatei Fourier discrete calculate în 512 puncte w=-pi:2*pi/512:pi-2*pi/512; plot(w,fftshift(abs(X))),grid %vezi help fftshift plot(w,fftshift(angle(X))),grid subplot(2,1,1), plot(w,fftshift(real(X))),grid subplot(2,1,1), plot(w,fftshift(imag(X))),grid Exerciţii Utilizând exemplul de mai sus, determinaţi şi reprezentaţi grafic partea reală, partea imaginară, modulul şi faza transformatelor Fourier discrete pentru următoarele semnale: = pentru 010n≤≤ [] [] 1xnnδ δ=−−−+− pentru 010n≤≤ [] [ ] [ ] [ ] 2135xn n n nδδ nπ ⎛⎞ pentru 020n≤≤ [] 3sinxn= ⎜⎟ 5⎝⎠ nπ ⎛⎞ pentru 020n≤≤ [] 4cosxn= ⎜⎟ 5⎝⎠ 6 3 Răspunsul la impuls al unui sistem liniar şi invariant în timp discret 6 3 1 Să se determine şi să se reprezinte grafic funcţia pondere a unui SLITD definit prin: 0 9 1 0 3 0 6 1 0 6 2yn yn xn xn xn−−= +−+− [] [ ][][][] Rezolvare b= ; a= ; [h,t]=impz(b,a) %s-a calculat răspunsul la impuls al sistemului Vezi help impz impz(b,a), grid stem(t,h), grid Verificaţi şi exprimaţi efectul comenzii: impz(b,a,40,3), grid Exerciţii Să se determine şi să se reprezinte grafic funcţia pondere a unui sistem definit prin: 1 [] 1 27 1 0 81 2 0 5 3 0 125 4yn xn xn xn xn xn=− −+ −− −+ − [] [ ] [ ][][] 0 9 1yn yn xn+−= 2 [] [ ] [] 0 13 1 0 52 2 0 3 3 0 16 0 48 1 0 48 2 0 16 3yn yn yn yn xn xn xn xn+−+−+−= −−+−−− 3 [] [ ] [ ][][][][][] 0 9 2 0 3 0 6 1 0 3 2yn yn xn xn xn+−= +−+− 4 [] [ ] [] [ ][] 66 4 Răspunsul unui SLITD la un semnal de intrare 1 Să se determine şi să se reprezinte grafic răspunsul unui sistem definit prin: 0 9 1 0 3 0 6 1 0 6 2yn yn xn xn xn−−= +−+− [] [ ][][][] =−− pentru 040n≤≤ la semnalul de intrare [] [] [ ] 110xn n nσσ Rezolvare b= a= x=[ones(1,10), zeros(1,31)]; y=filter(b,a,x); %vezi help filter n=0:40; subplot(3,1,1), stem(n,x), grid, title('x[n]') subplot(3,1,2), impz(b,a), grid, title('h[n]') subplot(3,1,3), stem(n,y), grid, title('y[n]') Exerciţii Să se determine răspunsul sistemelor definite prin relaţiile de mai jos 1 [] 1 27 1 0 81 2 0 5 3 0 125 4yn xn xn xn xn xn=− −+ −− −+ − [] [ ] [ ][][] 0 9 1yn yn xn+−= 2 [] [ ] [] 0 13 1 0 52 2 0 3 3 0 16 0 48 1 0 48 2 0 16 3yn yn yn yn xn xn xn xn+−+−+−= −−+−−− 3 [] [ ] [ ][][][][][] 0 9 2 0 3 0 6 1 0 3 2yn yn xn xn xn+−= +−+− 4 [] [ ] [] [ ][] la următoarele semnale de intrare: 1 [] = pentru 040n≤≤ [] 1xnnδ = pentru 040n≤≤ 2 [] [] 2xnnσ pentru 05n≤≤ 3 [] 3xnn= 7